REPORT No. 600 



AN ANALYSIS OF THE FACTORS THAT DETERMINE THE PERIODIC TWIST OF AN 
AUTOGEEtO ROTOR BLADE, WITH A COMPARISON OF 
PREDICTED AND MEASURED RESULTS 

By John B. Whbatlbt 



SUMMARY 

An analysis is presented of the factors thai determine 
the periodic twist of a rotor blade under the action of the 
air forces on it. The results of the analysis show that the 
Fourier coefficients of the twist are linear expressions 
involving only the tip-speed ratio, the pitch setting, the 
inflow coefficient, the pitching-moment coefficient of the 
blade airfoil section, and the physical characteristics of 
the rotor blade and machine. The validity of the analysis, 
was examined by using it to predict the twist of a rotor 
whose twist characteristics had previously been measured 
in flight. The agreement between the calculated and 
experimental results was satisfactory. An examination 
of the assumption used in the analysis — that the twist is a 
linear function of the radius — disclosed that the approx- 
imation introduced no appreciable error. From this ex- 
amination, a formula for the torsional rigidity of the 
rotor blade was derived. 

INTRODUCTION 

The development of the wingless direct-control auto- 
giro has heen hampered by a number of secondary 
difficulties. Probably the most troublesome are the 
avoidance of excessive or unstable center-of-pressure 
travel in the rotor and the elimination of rotor and 
control-stick vibrations. The production of a few- 
designs that are satisfactory in these respects has 
demonstrated that the difficulties are the designer's 
problem and are not inherent in the direct-control type 
of rotor; however, the large number of unsuccessful 
machines is evidence that the basic factors controlling 
the behavior of the rotor are as yet not clearly under- 
stood. 

A general survey of the problem indicated that both 
center-of-pressure travel and rotor vibrations are 
markedly affected by the periodic twist of the rotor 
blade arising from the interaction of air forces, elastic 
forces, and inertia forces during the flapping oscillation. 
It was accordingly decided that the factors controlling 
this twist must be understood before any real attack on 
the initial problems would be fruitful. This paper pre- 
sents an analysis of periodic blade twist in which the 



factors controlling the twist are studied. The analysis 
is supported by a comparison of predicted and measured 
twist on a direct-control type of autogiro. 

ANALYSIS 

The motion of an autogiro rotor blade consists chiefly 
of rotation about the rotor axis, oscillation about the 
flapping hinge, and oscillation in twist about the blade- 
span axis. Additional components of the motion are 
oscillation as a pendulum in the plane of the rotor disk 
about a second hinge and oscillation in bending in a 
plane containing the blade span and the rotor axis. 
Experimental evidence has shown that these additional 
components have only a second-order influence on the 
air forces acting on the rotor blade, and they will 
consequently be neglected in the subsequent discussion. 




Fiquee 1.— Force diagram of autogiro rotor-blade element. 



The coefficients of the air-twisting forces on an ele- 
ment of the autogiro rotor blade are diagrammed in 
figure 1. In general, the air forces on an airfoil will not 
pass through the aerodynamic center but will assume 
such a position that the moment of the air forces ex- 
pressed in coefficient form is constant about the aero- 
dynamic center. The component of the centrifugal 
force normal to the blade and the inertia forces of the 
blade pass through the center of gravity of the blade. 
Then the moment of the air forces about the center of 
gravity is the twisting moment on the blade. Let G u 
be the moment coefficient of the air forces about the 
center of gravity; then from figure 1 

(1) 
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where the symbols are defined by figure 1 ; Cr is positive when the center of gravity is farther from the leading 
edge than the aerodynamic center, and I is positive when the center of gravity is below the aerodynamic center. 

From the definition of moment about the aerodynamic center, if h is positive from the aerodynamic center 
toward the trailing edge, 

-O m c (2) 

where G„ is the coefficient of the moment about the aerodynamic center. Then 

Cu=O m -\~{Ci 1 cos ctr+Oo sin a r ) —\(C L sin a T —0 D cos a T ) (3) 

It will be assumed that a T is sufficiently small that the cosine differs negligibly from unity and the sine can be 
equated to the angle; also, G L will be assumed a linear function of the angle of attack. These assumptions do 
not accurately represent the conditions in the rotor when the resultant velocity is small and the angle of attack 
high, but the error introduced by them has been found to be small (reference 1). Then 

Cu=C„+a%{.a+C I >)-ka a *-C I >) (4) 

where a is the lift-curve slope (in radian measure). 

The lift-curve slope a at infinite aspect ratio lies between 5.8 and 5.9 for most airfoil sections. Inasmuch as 
0 D will have actual values varying from 0.009 to 0.03 below the stall, with a weighted average of about 0.015, it 
is thought unnecessary to use the G D term in the expression for G M . In addition, l/e will normally be less than 
0.02 and (a<z r a — G D ) will be of the same order of magnitude; the l/c term will consequently be dropped. 

From reference 1, if xR be substituted for r, the nondimensional velocity components at the blade element 
of an autogiro rotor traveling at a speed V equal to /iQB/cos a are: 

Ut=^£=x+h sin ^ (5) 



Up^^^X+^pXLx+^—fiac+xbi+^fio^ cos $+(— xai+~nb^ sin $+(j/Mi+2xb^ cos 2^ 
2xoj^ sin 2^-t~pOa cos Zf+^iib» sin Zf 



(6) 



where u T is the component of the resultant velocity perpendicular to the blade-span axis and to the rotor axis, u P 
is the component of the resultant velocity perpendicular to the blade-span axis and to Ur, and ip is the azimuth 
angle of the blade from its down-wind position. Also 

Ur s =x*+^fi 2 +2fi% sin i'—^ 2 cos 2ip (7) 

u T Up=x\+jv?b 2 +(^— vxa<> ^p&h^ cos a^— ^"J— \&b^ sin ^ 

+^^i+2a^6 4 ^ cos 2^+^—^fi 3 a 0 +iixb l —2x'a^ sin 2^+^— ^/rfy+^W) cos 3^ (8) 
+(j;Aii+f f<z&a) sin d^~^b 2 cos If+^Ot sin 4 ^ 

The acute angle /3 between the blade and the plane perpendicular to the rotor axis is described by the expression 

/3=Oo — <h cos ip — bi sin 4> — a% cos 2^— b<tsm2f — . . . (9) 

The notation of reference 1 will be used throughout this analysis; a list of symbols employed and of their 
definitions is given at the end of this section. 

The calculation of the air-twisting moment Mq at the hub end of the blade will be made on the assumption 
that the air forces he in a plane perpendicular to the blade span and depend only upon the resultant velocity in 
that plane. The angle <p between the resultant velocity and the plane perpendicular to the rotor axis will be 
assumed equal to its sine and tangent, and to have a cosine of unity. On this basis, the angle of attack of a blade 
element is 

*=*+*=g+* (10) 
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where 0 is the pitch angle of the blade, measured as the acute angle between the plane perpendicular to the rotor 
axis and the zero-lift line of the blade airfoil section. 

The pitch of the rotor blade will be the sum of the pitch setting and the instantaneous value of the angle of 
twist to which the blade is deflected by the twisting moment. If the pitch setting is given as 0 o +z0i and it is 
assumed that the twist is a linear function of the radius, the pitch angle 0 may be expressed in the form 

6=B 0 +xdi+xtQ+xe l cos f+xtji sin l/t+xet cos 2\f/+xq 3 sin 2$+. . . (11) 

The use of a Fourier series in ^ for the angle of twist is justified, as it was for the flapping angle f), by the 
fact that the twist angle must be a repeated function of ^. Some question concerning the assumption that the 
twist is linear along the radius naturally arises; the problem will subequently be examined in more detail. 

The air-twisting moment at the blade hub can now be expressed in integral form as 

M a = J'^p&BPuSCutHx (12) 

The integration is performed from B to 0 (where B is arbitrarily assigned the value l—e/2B) to allow for tip losses. 
Substituting for 0 M from (4) 

M Q =j o B ±pQWcw(c n +a* T ^dx (13) 

Expression (10) for a? can now be used 

Substitute for 6, Ur* and Ur u P ; then integrate and collect; and 

+ S[fe + ^ ! ' B ]) + ("^ MaoB3+6l [ l^+i^j-iMO^+^^+l^+f^-lM 2 ^) cos * 

-l^-^^+^^+^J+lwaB 3 -^^) cos 2* 

+(-^aoB-|-^i^-fa^+^^ sin 2^ 

+(-^ 1 5+|^-| / x s e 1 5 s -^'hB 3 +€{^*+^]) cos H 

+(|Ai5+|^ a B 2 -|Ai^+|^+'?{|^+^ J B 1 ]) sin 3*} 

In order to examine the variation of blade twist with radius, it is necessary first to establish an expression 
for the total air-twisting moment outboard of any station a;. This moment can be expressed simply as 



(16) 



(17) 
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Substituting as in (14) and integrating 

+[|6i+|M%][S 3 -a?]+| £ i[B 4 -a; 4 ]) cos I^JlB-*] 

[^-^+[-|^+|^+|w--gi«%][B«-afl +^ 1 [^- a; «])sin^+(M 8 [-^o-|^]{5- a; ] 
+[- JM 2 ^- i^] [IP-* 2 ] +[§&*- f MH.] t-S 3 -^ +j* [^-x 1 ]) cos 2^ 

+ (-^ 3 ao[5-a:]+[|Mii+^][5 a -^+[-|a 2 +^ £l - |^][B»- a ?]+i^[ J B»- a! »]) sin 2* 
+(-^ 2 6i [£-*]+ [1^- I^+I^^-^-I^^-^+l^f^-^cosS^ 

+(jM^[i^:e]+[§/«k-f sin 3* } 

The torsional deflection of the rotor blade can now be detecrmined. Let Q designate a twisting moment 
that is distributed along the radius in the same manner as Mq x and is of such magnitude that the blade tip is 
deflected through 1 radian; then denoting the instantaneous torsional deflection at the tip of the blade by v 

$^i+Gu=M Q (18) 

where g is the equivalent moment of inertia of the blade about the elastic axis and 

u=€o+«i cos ^+>?ismiH-^cos2#+>Jssin2^+' . . (10) 

The value of g may be arrived at by considering the problem to be that of the torisonal deflection of a bar 
connecting two masses having moments of inertia of infinity and zero. The system then has but one degree 
of freedom and only principal modes of vibration will be assumed existent. Rayleigh's approximation to tho 
exact solution of the problem (reference 2, p. 59), which is sufficiently accurate for use here, is 

■ 2=^ (20) 

where i is the moment of inertia of the bar per unit length, and I is the length of the bar. 

For a particular case, 1=0.0050, Z=20, <?= 1,700, and fl=21; then 2=0.0333 and fi s 2=14.7. It is seen 
that 4 Q*2 an< i even 9 Q 2 2 are quite small in comparison with O; consequently inertia effects on the torsional 
vibration will be neglected. Then 

(?{«o+«i cos iH-i7i sin ^+6j cos 2^+173 sin 2^+ej cos Zf+ri 3 sin 3^} 

4™^|#+|m s 3])+(-|^ cos i, 

cos 2^+(-|^aoB+iA^-3W+|^+, 2 [^+^]-|^) sin 2f 
+(-|m 2 6iB+|^B 3 -|aB 2 -^ 2 5 3 +^[|^+|^]) cos 3* 
+(|M a a l B+|^6 s B 2 -| /1 ^iB 2 +^^ 3 +773[|5*+|p 2 5 i ]) sin 3*) 



(21) 
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The expressions for the thrust and the flapping-motion coefficients have been established in reference 3 for 
the rotor with varying twist; they are: 

Cl = 2 ^ jx(g-l M »)+^|g+0a06,f)+^-^ (24) 

6l= ^^K^ +0 - 035M ^ a3B3 -5(l S5+ ^^)-i^} (25) 



(23) 



(26) 



It will be expedient to substitute average values for B and y in Oj and 6 a ; in actuality, 5 will be but little 
different from 0.970 for solidities near 0.05, and y will be between 10 and 18 for present rotor blades. The 
assumption that B is 0.970 and that y is 15.0 will accordingly introduce little error. The substitution will be 
made in such a way that the resultant expressions will be linear, and the coefficients of X, 6, e, and i\ will have the 
same form and exponents as the similar factors already present in the expressions for the twisting moment and 
the flapping motion. Examination of these equations discloses that the consequent forms for Os and b% are: 

a 2 =M 2 {0.08547\B+0.07467^oB 2 +0.0588T(ei+eo)B 3 } 



+ { i{-0.0184: ei +0.0Q12y Vl B i } + 0.0150ye i B !> +0.2200r} 2 B 

(29) 



b 3 = — M 2 { 0.306 \B" 3 + O.3820oB- a + 0.294 (0j + t^B~ l } 
[0.2200 1 „ „„™ „ , 2.650 



-e^+O.Ol^-O^OOe^+^^aB- 3 



Inspection of (26), (27), (28), and (29) shows that substitution for y 1 B i has been made in the denominator 
of (26) and (27), and in the numerator of (26) whenever y has been raised higher than the first power. The 
quantity in braces of (27) has been divided by B s and the resultant y^B* outside the braces has been evaluated. 

The solution of equation (21) for e and rj now will follow, after substituting for the a and b coefficients, by 
equating the coefficients of identical trigonometric terms in (21). Before this operation is performed, the work 
may be simplified to some extent by considering the order of accuracy required in the substitution. 

It has already been shown (reference 1) that the expressions for the thrust and torque are evaluated to a 
sufficiently high order of n (the fourth) if Oo, a*, and b 2 axe expressed to the order m s and a x and bi to the order 
H 3 . Reference 3 shows that the same order of accuracy for the thrust and torque will be obtained if ^ is evalu- 
ated to n*, 6i and iji to /x 3 , and ej and ij 3 to y?. The coefficients «3 and ij 3 do not, to the order of n 4 , influence the 
thrust and torque. Reference to equation (21) establishes that e n and tj„ are of the order p", a fact that has already 
been implied by the form of the expressions for dj and b 2 . 

It is seen now that Ci and &i may be expanded in a linear form which is developed only to the n z order and 
that all terms in cto above the order may be dropped. Then, substituting for (h and b 2 , 

a o =^XB 3 + 7 0 o (i^+| M ^)+ 7 (e 1 +e o )(^+iM 2 ^)+^7^ (30) 
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1.768 
— 0.853nej — ^grv"k 



& 1=M7 x(§iF-O.O542M 2 )+M70o(^+O.m 
+0.1676jt 2 7'7iB 3 +0.0100/x7«2£ 4 -0.853/n l2 



(3D 



(32) 



The term Mw/It^ 1 hi Oo has been neglected, since it amounts to 2 percent or less of Oo- 
Substitute for oo, ai, 6 1; a* and 6 2 in (21); there results a set of linear equations in e and v that can be solved 
in succession by starting with the highest order. Thus, from the coefficient of cos 3^, 



+ 0.01 Unye^E 1 — 0.1683/iIj-B 3 } 



(33) 



(35) 



(36) 



From sin 3^ 

^ 3 ^pQ^ca^0.27lg\+0.381p 3 W ( 34 ) 
From cos 2^, retaining only terms in y? or lower, 
<fc=ipfl 1 iPeac^0.796A+0 

From sin 2^, 

(?^=ipQmac r {-0.0291jtV>^-0.0290MV<>B s -0^^^^ 
— O.OlOOT^+O.lOSSijiB 4 } 
From cos ^, retaining all terms up to /x 4 , 

- £a (iB*-0.073M 2 5 2 )-0.0069/x J T'7i5 9 -0.0005MTe2B 7 -0.0057/iJj 1 B 3 J 
From sin ^, 

+0-006 ^B3_0g3 W2+ ^} 
• From the constant term 

+ipnwt7 m (|s 3 +|^) ■ 



(37) 



(38) 



(39) 
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Comparison with (22) shows that the first part of (39) differs only by insignificant terms from the thrust 
multiplied by cct/o-ttR, or, since <r is bc/nR, 



(40) 



Let pcatfR^CrlZQ be put equal to A; then after substituting from (35), (36), (37), and (38) for ej, r^, e,, and 17, 
and neglecting insignificant terms, the final expressions for the twist coefficients are obtained: 



£l = - w a(x(^£*-O.O161m s If)+0o(^^ 

+^74 s (o.0021XB 7 +0.00075 tt B 8 -0.0002(d 1 +^)B 9 +0.0071^B 8 ) 
I ac T J 

1Jl =M{xQ-B+O.34l£)+fl o ^5 i +O.233^)+O.175M s (0 1 + e o)5+^B 2 J 



003 



4V^ 2 {o.OO7\B 3 +O.OO60 o B 4 +O.OO5(0 1 +<*).B 5 -O. 

e i =^A\0.7Q6\+0.57Se o B+0.55Ue 1 +e o )B i -l^B\ 
I 2 OCt J 

~u. s A 1 \om2\B i +Q.OUe a B l +Om9(6 1 +e t )B< s -0.055^B i 
I ac T 

172=— Ai 2 Y-4{0.0291XB 4 +0.0290e () B 5 +0.0225(e 1 +e 0 )5 fl } 



-fJyA* 



e 3 =M 3 7^{o.0085XB 3 +0.0143e 0 B 1 +0.0108(fl 1 +«b)B 5 J 

+M 3 TA 2 fo.O135X J B 7 +O.Ollie o B 8 +O.OO97(0 1 +e o )5 9 -O.OO47^5 8 l 
1 acr J 



0.271, 



fls=M , 4(-^pX+O.381flb+O.28O(0i+$O5 

+^A 2 (o.098\B 3 +0.048eoB*+0.049(e 1 +eo)5 5 -0.023^5 4 
I OCt 



The expressions for the twist coefficients disclose that 
only <=o, rj u and ^ involve the factor A^^> which reduces 

to ,jg pcWRPOn and is independent of c T . The rest 

of the coefficients — ei, 170, «3, and 173, and parts of «o, 
771, and 63 — are proportional to A and consequently to 
c T . Thus if c T is zero, only «o, iji, and ej differ from zero, 
and then only if O m is not zero. Exclusive of the 
moment arising from O m , the factor A represents in 
nondimensional form the ratio of the moment of the 
air forces to the torsional rigidity of the rotor blade. 

The probable magnitude and range of values of A 
can be estimated. The rigidity O will be proportional 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



to the polar moment of inertia of the blade spar and 
inversely proportional to the length of the blade. The 
moment of inertia of the spar will be proportional to the 
fourth power of the blade thickness. Then 



0.0078 \B*+ O.OO570 O B 9 +O.OO52 (», + toW°- O.OOSO^B 9 ) 



OCt J 



ffoc 



~R~ 



(48) 



where t is the blade thickness divided by the chord. 

The numerator of A can be examined by the following 
considerations. In the design of a rotor the pitch 
setting chosen is almost invariably the one that results 
in the highest efficiency. The rotor speed is then 
adapted to varying maximum speeds by adjusting the 
solidity. The rotor disk loading is fixed between fairly 
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narrow boundaries by the requirement of good low- 
speed performance. For a given pitch setting, C T /<r 
(ratio of thrust coefficient to solidity) is constant and 
C T pQ?R 2 is equal to the disk loading. Thus, aptfR 2 -will 
be almost constant. Now, from (48) and this discus- 
sion, 



m* 



Or, from c=- 



7tR 



Acc 



EM 



(49) 



(50) 



It has been found that a value of A of 0.600 is 
associated with a rotor that has a radius of 20 feet, a 
chord of 1.00 foot, a blade thickness ratio of 0.175, and 
c T equal to 0.038 foot. It seems unlikely that R/d will 
increase by more than 25 percent above the value given 
here of 114; assuming this increase, and an increase in 
Cr/c from 0.038 to an upper limit of 0.06, A would 
become approximately 1.90. The lower limit is ob- 
viously zero, since c T jc may become zero. It will be 
found that normal designs will result in a value of A 
of less than unity and that the given value of 0.600 is 
larger than the average. 

LIST OF SYMBOLS 

R, blade radius. 

b, number of blades. 

c, blade chord. 

c T , distance between aerodynamic center and 

center of gravity of rotor-blade element, 
r, radius of blade element. 
x, r/R. 

6 0 , blade pitch angle at hub, radians. 
0 1; difference between hub and tip pitch angles, 
radians. 

coefficient of cos ruff in expression for 0, 
radians. 

•tin, coefficient of sin nty in expression for 0, 
radians. 

0, instantaneous pitch angle, radians. 

4>, blade azimuth angle measured from down 

wind in direction of rotation, radians. 
Q, rotor angular velocity, d*j//dt, radians per 

second. 

speed of axial flow through rotor. 
liQR, component of forward speed in plane of disk, 
equal to V cos a, where V is forward speed, 
feet per second. 
0, blade flapping angle, radians. 
da, coefficient of cos Tup in expression for j8, 
radians. 

b n , coefficient of sin nip in expression for 0, 
radians. 

I u mass moment of inertia of rotor blade about 
rotor hinge. 



a, rotor angle of attack, radians. 
C m , pitching-moment coefficient of rotor-blade 

airfoil section. 
M r , thrust moment about horizontal hinge. 
M w , weight moment of blade about horizontal 

hinge. 

Ui&R, velocity component at blade element per- 
pendicular to blade span and parallel to 
rotor disk. 

-Up&B, velocity component at blade element per- 
pendicular to blade span and to u^iR. 
T, rotor thrust. 
C t =T/pQ s tR* 
a, slope of curve of lift coefficient against angle 
of attack of blade airfoil section, in radian 
measure. 

¥>=tan -1 — 

a T , blade-element angle of attack, radians. 
y= j ?mass constant. of rotor blade. 

_B=1— 2"g> factor allowing for tip losses. 

67, torsional rigidity of rotor blade, ft.-lb. per 
radian. 

FLIGHT TESTS AND CALCULATIONS 

Data for investigating the validity of the analysis 
were obtained in flight tests of a Kellett KD-1 direct- 
control wingless autogiro. The physical characteristics 
of the machine and its rotor are given in table I. Meas- 
urements were made in a steady glide of the air speed, 
the rotor speed, and the blade motion. The air speed 
was obtained with a trailing pitot-static head and an 
N. A. C. A. air-speed recorder; the rotor speed was ob 
served with a calibrated rotoscope; and the blade motion 
was photographed withamotion-picturecameramounted 
on and toning with the rotor hub. The photographs 
obtained with the motion-picture camera established 
the blade flapping motion and, in addition, the instan- 
taneous twist at the rotor radius of the markers on the 
blade. 

TABLE I. — PHYSICAL CHARACTERISTICS OF KD-1 
AUTOGIRO 

Gross -weight, W 2,100 pounds. 

Rotor radius, R 20.0 feet. 

Blade weight, v>t 61.5 pounds. 

Blade-weight moment, M w 482 pound-feet. 

Blade moment of inertia, Ii 176 slug-feet*. 

Blade chord, c 1.00 foot. 

Chordwise location of blade center of gravity 

from leading edge 0.280 foot. 

Chordwise location of aerodynamic center 

from leading edge 0.242 foot. 

Number of blades, 6 3. 

Rotor solidity, <r 0.0478. 

Blade airfoil section GSttlngen 606. 

Blade pitch setting (constant), 0 0 0.0960 radian. 

Airfoil section moment coefficient, C„ (about 

aerodynamic center) —0.056. 

Blade torsional rigidity constant, 0 1,700 pound-feet. 
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In order to investigate the validity of the analysis, 
the flight-test data were used in two ways. The analy- 
sis was checked directly by predicting the blade twist at 
0.75 B from the physical constants of the rotor and the 
value of the inflow coefficient X. The factor X was 
calculated from the experimental thrust coefficient by 
the substitution of known values in the expression for 
the thrust coefficient given in the analysis. A further 
examination of the analysis was made by substituting 
the experimental values of the inflow factor, the blade- 
motion coefficients, and the twist coefficients at one tip- 
speed ratio in the equation expressing the twisting 
moment as a function of the radius. The resultant 
twist deflection for a blade of constant rigidity followed 
directly and could be qualitatively compared with the 
basic assumption of a linear variation of twist with 
radius. 

RESULTS AND DISCUSSION 

Measured values of rotor speed and thrust coefficient 
are shown in figure 2, and derived values of the inflow 
factor X are given in table II. The experimental blade- 
motion coefficients are presented in figure 3. 
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Fioube 2.— Hotor speed and throst coefficient of KD-l autoglro rotor as measured In 
flight; p= 0.00231 slug/cu. ft, 

TABLE n.— DERIVED VALUES OF INFLOW FACTOR X 





X 


f 






X 


a 125 

.160 
.176 


0.0214 
.0213 
.0212 


0.200 
.226 
.260 


0.0209 
.0204 
.0197 


a 276 
.300 
.326 


a 0191 
.0186 
.0179 



The twist coefficients shown in figure 4 represent a 
comparison of the experimental points with the cal- 
culated values for the radius at which the measurements 
were made. The agreement of theory with experiment 
is satisfactory, and strikingly so for the largest 
coefficient. The calculated eo is consistently smaller 



than the measured value, but a reasonable explanation 
of the disagreement will be given later. Unfortunately, 
the experimental results for ti and ea are badly dispersed; 
the mean of the points is not appreciably different, 
however, from the predicted values. 

The variation of twisting moment and twist angle 
with radius is illustrated in figures 5 and 6. The 
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Fionas 3.— Flapping motion coefficients of KD-I autoglro rotor as measured In 

flight. 

twisting moment at the hub for each component of the 
twist has been considered unity; the blade rigidity 
from a;=0 to £=0.05 has been considered infinite as 
an approximation to the high rigidity inboard of the 
vertical pin; and the rigidity outboard of 3=0.05 was 
chosen to make eo unity at 3=1.00. Examination of 
the figures discloses that eo attains a larger value at 
3=0.75 for the same twisting moment at the hub 
than any of the remaining components of twist. This 
result can be considered a partial explanation of the 
underestimation of eo in figure 4. The curves in figures 
5 and 6 further indicate that the assumption of linear 
twist is a reasonably accurate approximation to the 
actual variation. 

The curves in figures 5 and 6 suggest that it would 
be erroneous to calculate the torsional rigidity from 
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FiotntE 4.— Measured and calculated twist coefficients of KD-1 antogiro rotor. 
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the twist produced by a constant moment along the 
blade. Instead, it is recommended that a moment be 
assumed which has a value Q at the blade root and 
varies with the radius according to the expression 
(B 1 — X s ). The definition of 0 suggested is that it be 
the moment at the hub, distributed as indicated, which 
will produce a twist of 0.80 radian at 0.75 B. Thus, if 
I„ is the polar moment of inertia of the blade cross 
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Fionas 6.— Eadlal distrlbntlon of twisting moment and twist angle of KD-1 autoglro 
rotor; sin mi> components; p—0J321. 

section, B the radius (blade length), and E, the modulus 
of elasticity in shear, 

„ 64 EJ, 
585 B 

where O is in pound-feet. 

E, is in pounds per square inch. 
Ij, is in inches*. 
B is in inches. 
The value of 0.80 rather than 0.75 radian for the deter- 
mination of O appears to result in a curve that is 
better approximated by a straight-line distribution to 
1 radian at the tip, as evidenced in figures 6 and 6. 

The merit of the analysis made in this paper depends 
upon the accuracy and facility with which it may be 
used to predict the twist of a rotor blade before the 
rotor itself has passed the drawing-board stage of 
design. The use of the analysis in this manner is not 
obvious since at first it appears that the rotor speed is 
required in order to calculate the twist; whereas the 
calculation of the rotor speeds can be made only after 
the twist, and consequently the thrust coefficient, is 
known. In one sense this objection is valid but in 
another the difficulty mentioned is not insurmountable. 
Assume a rotor design in which the known factors are 
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the blade moment coefficient, the distance between the 
aerodynamic center and the center of gravity, the radius, 
the chord, the torsional rigidity, and the design maxi- 
mum speed. Because of the variation of rotor effi- 
ciency with tip-speed ratio, it is mandatory that at 
maximum speed the tip-speed ratio shall be between 
0.40 and 0.45. When the tip-speed ratio for design 
maximum speed is chosen, the tip speed is fixed. The 
twist coefficients e 0 , vu and which depend principally 
upon O mi can now be found with satisfactory accuracy 
by using the values of pitch setting and inflow coeffi- 
cient X that would be assigned to the rotor if the twist 
were zero. For a given airfoil section there is a mean 
lift coefficient which results-in maximum efficiency and 
which fixes the ratio 0 T /a; tho values of X and the pitch 
setting corresponding to this mean lift coefficient should 
be used. The coefficients e 0 , j? u and cj now are known 
and fixed, and their effect upon C T fo- can be evaluated. 
The desired value of 0 T /<r is now attained by adjusting 
the pitch setting to offset the effect of twist, and all 
design requirements have been met. A final check of 
the twist coefficients using the final values of X .and the 
pitch setting can be made but, since O m dominates the 
only coefficients that influence C T l<r, it is found that only 
in exceptional and peculiar designs will the twist 
coefficients be affected. 



CONCLUSIONS 

1. The assumption that the twist of a rotor blade 
varies linearly with the radius is a satisfactory approxi- 
mation to actual conditions. 

2. The analysis of. blade twist predicted without 
important error the twist of a rotor used as an example. 

3. The torsional rigidity of the rotor blade should 
be calculated on the basis that the twisting moment 
varies with radius as (B^—x 1 ). 
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